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ABSTRACT 


This report examines the possible use of "rough,” or nondiffraction- 
limited, reflectors for collecting optical signals. It is shown that in the 
absence of background radiation, the reflector's surface quality has little 
effect on the performance of a properly designed receiver, but that the pre- 
sence of even small amounts of background radiation can lead to severe perfor- 
mance degradation. Techniques are suggested for improving receiver performance 
in high-background environments, and bounds and approximations to the exact 
error-probability expressions are derived. 
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SECTION 1 


INTRODUCTION 


The development and eventual deployment of large, Earth-orbiting 
optical receivers is under consideration at the present time. Such receivers 
could be used to detect weak optical signals from distant deep-space vehicles, 
or to relay high-rate optical data from nearby planetary probes (and other 
near-Earth sources) to the ground. From a communications viewpoint, one of 
the most important parameters of an optical receiver is its effective aperture, 
which determines the amount of signal power the receiver can extract from an 
incident signal field. Dif fraction- limited optical reflectors provide the 
greatest protection against background interference, allowing the minimization 
of the receiver's f ield-of-view . However, it is generally both difficult and 
costly to construct large, diffraction-limited reflectors at optical wave- 
lengths, hence our interest in the potential use of "rough" reflectors for 
optical communications. Thus, multielement reflectors designed for coherent 
operation at infrared wavelengths could be employed in a "photon-bucket" mode 
at much higher frequencies, or large, inexpensive reflectors with poor optical 
surface quality may at times be employed to collect signal energy. Since 
"rough" reflectors typically scatter a significant fraction of the incident 
signal field, special processing techniques must be used to recover the 
scattered signal energy to improve receiver performance. 

The following sections address some of the fundamental issues 
associated with the use of nondiffraction-limited reflectors for optical 
communications. In Section 2, a useful model for signal and background 
intensity distributions is developed. The structure and performance of a 
minimum error-probability, direct-detection optical receiver is derived in 
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Section 3, along with bounds and approximations to the exact expressions. 
Finally, in Section 4, numerical examples are provided to illustrate receiver 
performance under various operating conditions of interest, and to provide in- 
sight into the use of bounds and approximations developed in the previous 
section. 
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SECTION 2 


OPTICAL RECEIVER CHARACTERISTICS 


The class of optical receivers considered in this report consists 
of the following basic optical components: a reflector (of radius R and focal 

length f) designed to collect and focus optical fields, a narrowband optical 
filter (passband [A, A + AA)) designed to pass the signal fields but 
reject out-of-band background radiation, a spatial filter to limit the 
receiver's f ield-of-view , an optical detector array to convert the impinging 
optical fields into electrical signals, and a post-detection processor to 
attempt to reconstruct the transmitted message with the desired fidelity. The 
arrangement of these components is shown in Figure 2-1. 

The reflector's surface can be modeled as an isotropic random 
function with normal height distribution above and below the mean surface. 

Let the random variable h denote the deviation of the actual surface from its 
mean value at a particular point on the reflector, and let be the 
distance between two points on the surface. It is generally assumed [Refs. 1 
and 2] that h is Gaussian with probability density 

p(h) = — i— e~ h /2 ° (2.1) 

/ 2tt a 

and that the correlation function along the surface is of the form 


C ( T ) 


C 




e 


( 2 . 2 ) 
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Figure 2-1. Optical Receiver Components 



Thus, a is the root-mean-square (rms) surface deviation, and T is the 

c 

correlation length averaged over all directions along the surface. 

A point in the focal plane can be specified by its Cartesian 
coordinates (x,y), or, equivalently, by its polar coordinates (0,<|>) (for small 

r~2 2 a a 

0,0 - vx + y /f and (j) is the angle between the vector _r = xi + yj and one of 
the Cartesian axes). If a circular reflector with rms surface deviation a and 
correlation length is illuminated by a normally incident plane-wave, then 
the fraction of power passing through a circular aperture of radius 0f 
centered in the focal-plane is given by the expression [Ref. 2] 


f(0) = f L ( 0 ) + f 2 ( 0 ) 


(2.3a) 


where 


F^O) = e g 


jl - ^J o 2 (^R6) + J^^Re)^ | , 


(2.3b) 


F ( 0 ) = 1 - e 


1 + 


-G*T c e/2) 2 /m 

I ^ 


m: 


m=0 


( 4tt a \ 

* ) 


(2.3c) 


(2.3d) 


(A is the center wavelength of the optical signal, and 6 A 2tt/A.) It is apparent 
that F^(0) is the contribution of the Airy pattern generated by the mean 
reflector surface, attenuated by an exponential function that depends on a, 
while F^(0) describes the distribution of power scattered out of the Airy 
pattern (note that F^(0) is independent of the reflector radius). For a 
"rough reflector" g >> 1 , and F(0) can be approximated as [Ref. 2] 


F(0) 


1 - e 


i ( 6 /e o > 2 


; g » i 


(2.4) 


where 
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0 n A /2 o/T 
0 — c 


The behavior of F(6) is examined in Figure 2-2 for various combi- 

nations of R, o, and at a wavelength of A = 1 pm. For g << 1 , the enclosed 

power depends strongly on the reflector radius R, as shown in Figure 2-2(a). 

As a (and hence g) increases, the contribution of the Airy pattern becomes 

less important (see Figure 2— 2( b ) ) * For large values of g (corresponding to 

a = 0.2 ym) , the contribution of the Airy pattern becomes negligible and the 

approximation of Equation 2.4 may be used. Note that for a given rms surface 

deviation, the amount of power scattered out of the Airy pattern does not vary 

significantly as the correlation length changes, but the actual distribution of 

scattered energy does (see Figure 2-2(c)). 

Equation 2.3a can be used to relate the power reaching the detector 

array to the intensity of the optical field at the receiver. If the intensity 

2 

of a normally incident plane-wave at the receiver is W/m , then the power 

passing through a circular aperture of radius 0f in the focal plane is 

P (0) = A 1 F( 0) (2.5) 

s r r 

2 

in watts, where A At tt R is the effective aperture of the reflector 
r — s 

(t is the bulk power transmission factor of the receiver, with 0 < t < 1). 
s s 

Background radiation may also enter the receiver along with the desired signal 

field. This type of radiation is generally of thermal origin, and it may be 

due entirely to self-emission by hot sources (such as stars) or it may be a 

combination of self -emission and reflected radiation (as from planets). With 

the exception of the Sun, stars generally appear as point sources due to their 

great distance from the receiver, while planets appear as distributed sources 

to receivers capable of resolving their disks. Here we consider only extended 

background sources that are characterized by a spatially constant spectral 

2 

radiance function, N( A) (the units of N( A) are W/m -ym-sr). The power 
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Figure 2-2. Fractional Signal Power F(0) as a Function of 0 for 

Various System Parameters: (a) a = 0.03 yin, T c = 0.25 m, 

and X = 1 ym 




Figure 2-2. Fractional Signal Power F(0) as a Function of 0 for Various 
System Parameters: (b) R = 5 m, T c = 0.5 m, and X = 1 ym 
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Figure 2-2. Fractional Signal Power F(0) as a Function of 0 for Various 
System Parameters: (c) a = 0.07 ym, R = 5 m, and X = 1 pm 






distribution in the focal plane due to a spatially constant background source 
is also constant; hence the background power passing through a circular 
aperture of radius 0f in the focal plane (over the wavelength range [X, X + AX)) 
becomes 


p b (0) 




TT 0 


0 « 1 rad 


( 2 . 6 ) 


in watts. Note that the background power reaching the detector array continues 

to increase with 6 , whereas the power due to a normally incident plane-wave 

reaches its maximum value as F(0) approaches one. 

An optical detector suitable for photon-counting applications 

generates an electrical signal that can be related to the number of photons 

detected in a given time interval. If an optical plane-wave of intensity 1^ 

2 

W/m is normally incident on the reflector, then the average number of 
photons detected by a centered , circular detector element of radius 0f over 
the time interval [t, t + At) can be expressed as 


K (0) = I F( 0) 
s s 


(2.7) 


with 


I 

s 




At 

r 


where I is the average number of signal photons passing through the 
receiver's effective aperture in At seconds, n is the detector quantum 
efficiency, and hv is the energy of a single photon. Similarly, the average 
count generated by multimode background radiation over the same time interval 
in the wavelength range [X, X + AX) becomes 


2-8 



( 2 . 8 ) 


K b (6) - i b . e 2 


with 


: b = 


-r— A At 
hv r 


/ 


X-f AX 

N(c)dC 


in counts/ steradian (sr) . As an immediate extension of the above results, note 
that the average count generated by a centered detector-ring of outer radius 
and inner radius Q^f is simply K( - K(0^)« This expression holds for multi- 
mode background fields as well as for normally incident plane-waves. In the 
following sections, the above models will be applied to determine the perfor- 
mance of optical receivers employing rough reflectors, and to develop tech- 
niques for improving receiver performance in the presence of interfering back- 
ground fields. 
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SECTION 3 


OPTICAL RECEIVER PERFORMANCE 
3.1 DECODER STRUCTURE 

In order to transmit information over the channel, the transmitter 
modulates the intensity of the transmitted field, which results in a corre- 
sponding modulation in the intensity of the received field. For long-distance 
communications, the received field can be modeled as an intensity-modulated 
planewave. The receiver bases its decision on a fundamental array of observ- 
ables , called counts , generated by each element of the detector array over a 
suitable set of time intervals, in response to the modulated received field. 
Each time interval should be chosen small enough to ensure that the intensity 
of the received field remains essentially constant over that interval. Con- 
sider a situation in which the transmitter sends one of M waveforms during 
every signal interval. Then, if each of M received symbols lasts exactly T 
seconds , N photodetectors responding to both signal and background fields over 
L time intervals (whose union comprises the interval [0,T)) generate the count 
array 



fk 

**• k lL 

[k] = 

k 

* * ' k NL_ 


(3.1) 


where, under hypothesis "i" (or H. ) , the count k ft is a Poisson distributed 
y r l in 

random variable with average value 


K 


£n 


= K 


1 ,s 
in 


+ K 


£n 


(3.2) 


3-1 



is b 

(Here K * refers to the average count due to the signal under H., and K is 
£n & l in 

t* Vi 

the average count due to the background radiation over the n detector element 
t h 

and £ time interval. The dependence on 0 has been suppressed for the sake 
of notational simplicity.) 

Since Poisson counts generated by a deterministic signal field or 
multimode background radiation over disjoint space (or time) intervals are 
independent random variables, it follows that the joint probability density of 
the elemental counts, conditioned on hypothesis "i" being true, can be 
expressed as 

k„ 


P([k] | H.) = exp 


L N \ L N /„i 

i \ T~r 1 r \ £n 


2_,2>w n n 


£n 


£=1 n=l 


£=1 n=l 


k ! 
£n 


(3.3) 


The M maximum a-posteriori M (MAP) decoder selects the hypothesis with the 
greatest probability of occurrence, given the observable array [k] . Using 
Bayes" rule and keeping only terms that depend on , the decisions can be 
based on the set of log-likelihood functions , where each element of 

this set is defined as 


D. A 

l — 



(3.4) 


Here the truncator keeps only hypothesis-dependent terms, 

observable array, the MAP decoder computes (fh} and selects H 
D = max {D.} 

q , i 


For each 
if 


L N 


max ^log^P(H^) 
i 


K 1 + 
s 



0 




l=\ n=l 


where i = 1, 2, M. 


(3.5) 
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In case several decision functions are tied for largest, the 

decoder may employ any predetermined strategy to choose one of the candidate 

hypotheses. In Equation 3.5, K. 1 and are defined as 

s £n 


L N 



( 3.6a) 


(3.6b) 


Therefore, the MAP decision is based on the a priori probabilities P(1L) , the 
total signal energy , and a weighted sum of counts obtained from every 
detector element over the entire observation interval. 

The decision function simplifies for the case of equilikely, equal- 
energy M-ary pulse-position modulated (PPM) signals. For this modulation for- 
mat, a T-second optical pulse in the i ^ T-second time slot (t _A T/M) 

represents hypothesis "i" [Ref. 3]. Now L = M, P(H.) = 1/M, and K 1 = K (since 

1 s s 

is s 

K&n = » independent of "i")* Furthermore, it can be seen from Equations 

3.2 and 3.6b that if does not depend on ,, £ n , then 

£n ’ 


i 


w 


£n 



i=£ 


i^£ 


Hence the decision functions reduce to 


(3.7) 
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(3.8) 


N 

D. = y k. w 
1 in n 

n=l 

where i = 1 , 2 , . . . , M. 

In this case, the decision is based on the sum of weighted Poisson 
counts from all detector elements over each time interval [(i-l)x, ii) , and the 
hypothesis corresponding to the time interval containing the greatest weighted 
sum is selected. 

Since each decision function is a discrete sum, several decision 
functions may take on the same maximum value with nonzero probability. When 
this occurs, the problem can be resolved by making a random choice among the 
decision functions, or by any other strategy that selects a hypothesis 
corresponding to one of the largest values. Decoder performance is not 
affected by the particular strategy employed to resolve ties. 

3.2 DECODER PERFORMANCE 

Prior to observation, each decision function is a random variable, 
since the individual counts are random variables. The statistics of the counts 
depend on which hypothesis is actually true. Given that is true, the 
probability of a correct decision is the probability that D^ exceeds all other 
D^ (i ^ q), plus the probability that the correct hypothesis is selected when 
several decision functions are tied for largest. In order to compute the 
average probability of a correct decision, the probability density of each 
decision function must be known. If the conditions of Equation 3.8 are 
satisfied, then the conditional probability density of the i^ decision 
function, D^ , given that H is true, can be expressed as 
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P D. (a lY = 2-J P k (i,H q ) 6(a ' V 
1 k=0 


(3.9) 


As explained in Appendix A, p (i,H ) is the probability mass over the point 

x q 

a^, when is true and when the weights are according to H^. 

When N detectors are used to observe M-ary symbols, the average 
probability of a correct decision is 


N P M (C> - 


E 


nV c 'V 


(3.10) 


If equilikely, equal-energy M-ary PPM symbols are transmitted, then by symmetry 


» P M (C) ' N P H <C 'V 


E fey (v) E 


P k (<,.H q ) p R <i,H q > 


1 E 

-l _n 


Pj (i-v 


M-r-1 


+ M (P 0 (q,H q ) 


i*q / 


(3.11) 


If a single detector is employed, then for Poisson distributed 


counts. Equation 3.11 becomes 


l p M (C) - 


M-l » 

E fey n E 


(K s+ K b ) k -(K s+ K s ) K b k -K b ] r 
k! e — e 


K j 

K b 


E 


V 


(M-r-l) 


' ( V M V 

+ M e 


(3.12) 


(K and K are the average signal and background generated counts, respectively, 
S D 

in each x-second time interval.) Note that as approaches zero, the average 


symbol error probability reduces to 
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lV E > ■ 1 - 1 P M (C > - (¥) 


e 


(3.13) 


-K 

s * s * 0 

which is just the probability of erring when attempting to decide between M 
empty slots by random choice. 


3.3 BOUNDS AND APPROXIMATIONS 

It is apparent from Equation 3.2 that the average count generated by 
any detector element is always at least as great as the average count induced 

by the background radiation. As the smallest element of [K^ ] increases, the 

£n 

envelope of the density of each elemental count approaches a Gaussian density 
with mean and variance both equal to the particular count-average [Ref. 4]. 
Thus, for large background levels, each detector output can be modeled as a 
Gaussian random variable. Since the sum of independent weighted Gaussian 
random variables is also Gaussian, the envelope of the probability density of 
the weighted sum defined in Equation 3.8, approximated as a continuous 
Gaussian random variable, becomes 

1 2 2 

? - 7 -(a - a. ) / 2af 

P_ (a|H ) * ( 2 tt of) e 1 1 (3.14) 

D. 1 q l 

l 

where, for equilikely, equal-energy M-ary PPM signals 


a . 

l 




i * q 


i = q 


(3.15a) 
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and 



(3.15b) 


b s 

In the above equations, K is the average count due to background radiation, K 

n n 

the average count due to signal , and the weighting factor associated with 
the detector element. Ignoring count equalities, the average probability 

of a correct decision becomes 


N P M (C) " N P M< C IV * 


f 


- ( a- a q )^/2 a q 2 J />a -(g-a.) 2 /? a. 2 

da — < / dp' 


(M-l) 


/ 2tt a 2 

q 


/ 2tt a 2 
i 

i * q 


(3.16) 


which yields, after some manipulation, 


n p m ( e > 


/ 


i-l dy 


rv 2 A 


v/2^ 


(M-l) 


a ( a -a. )) 

q „ , q 1 


1 - — Erfc I — L y + 

2 la, ; a. 

1 1 

i / q 


(3.17) 


This last expression is considerably easier to evaluate numerically than the 
exact expression given by Equation 3.11, and can therefore be used to approxi- 
mate the average symbol-error probability when the intensity of the background 
radiation is sufficiently great. The accuracy of this approximation under 
various operating conditions will be examined in the following section. 
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The performance of the single-detector receiver observing 


equilikely , equal-energy M-ary PPM signals can be bounded by the inequality 


P (E) < (M-l) P (E) 
1 M 12 


(3.18) 


which is proven in Appendix B< 


The single-detector binary error probability can be expressed as 


[Ref • 5' 


i p 2 <E > ■ e 


£ 


F n ( *> - i V” 


(3.19) 


where A = / K + K, - /IT , i|» = 2/(K + K,)K , and F (i|>) = e (<|>). Since 

sdd sbb n n 


for any ty > 0 and n > Q> I^C^) * I ( ^ ) > it follows that 


1 P 2 (E) < e 


oo 

£ 




li V*> 


- a 2 \A . 1 


= e “ — t— • + -y F (♦) 

A 2 o 


Al P 2 


The upper bound on the binary PPM error probability, ^P^ , can itself 
be approximated accurately in the limit of small or large K^. Thus, for 
K b « 1, 

i p 2’K i2 ; s«* (3 - 2u: 


while for K, >> 1 , we have 
b 




/ 2tt ^ 


• K b » * 
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The approximation in Equation 3,21b generally yields accurate results even for 

1). These approximations can be easily 
evaluated with a calculator to provide the designer with quick, accurate 
estimates of receiver performance. 


moderate values of (that is, 
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SECTION 4 


NUMERICAL RESULTS 

4.1 RECEIVER PERFORMANCE 

In Section 3 we showed that the error probability expressions (as 
well as the bounds and approximations) depend on the average signal and back- 
ground counts generated by the detector elements over the signal interval. For 
a given wavelength (henceforth we assume that X = 1 pm), the average signal 
and background counts, K g and , in turn depend on the reflector parameters 
a, T , and R, and on the angle 0 (which defines the receiver's f ield-of-view) , 
as was shown in Section 2. In general, 0 has to be optimized to achieve best 
performance (that is, lowest error probability) as the various system param- 
eters are changed. Next we examine the performance of an optical receiver 
observing M-ary PPM signals, for various reflector parameters and background 
environments . 

4.1.1 No Background Radiation 

First consider the limiting case = 0 , and define K gQ as the 

average signal count over the entire detector-plane. In this case no advantage 

can be gained by using multiple detectors; hence, we need only consider the 

single detector case (N = 1). Now ^P^(E) is a function of R, a, T , M, 

and 0. The variation of kT P k ,(E) with 0 for several values of a is shown in 

N M 

Figure 4-1, with R = 20 m, T = 1 m, M = 32, and K = 10. The limiting error 

c so 

probability of ( (M - 1 )/M)exp(-K go ) can be approached arbitrarily closely by 
increasing 0 for any choice of a; however, it is approached most rapidly by 
the diffraction-limited reflector (a = 0.0pm). Therefore, in the absence of 
background radiation, the only disadvantage of poor reflector surface quality 
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is that the larger f ield-of-view required to achieve a given performance level 
increases the chance that nearby stars, or other background sources, might 
enter the receiver's f ield-of-view. 

4.1.2 The Effects of Background Radiation 

As the background level increases from zero, optimization over 0 

becomes necessary to minimize the error probability. (The units of are 
-12 

counts per 10 steradians, or counts/pico-sr , referenced to an arbitrary 

timeslot. Hence, the average noise-count observed over any M-ary time-slot is 

2 

taken to be K, = I tt 0 .) This relationship between 0 and I is illustrated 
b b b 

in Figure 4-2, where the optimum 0, 0 , is displayed for various background 
levels, with the other relevant parameters held fixed (K gQ = 10, R = 20 m, 

T^ = 1 m, a = 0.1 pm, and M = 32). The dependence of 0^ on is quite evident 
from this figure: when the background level is increased, the optimum angle, 

0 q , decreases as the receiver attempts to exclude the interfering background 
radiation while collecting signal energy. However, the minimum error probabil- 
ity also increases with background level, because even with an optimized f ield- 
of-view, less signal energy is collected in a noisier environment. 

The minimum symbol error probability is shown in Figure 4-3 as a 
function of K gQ for various background levels and rms surface deviations 

(R = 20 m and T = 1 m) . A single circular detector of radius 0 f is used in 
c o 

both cases (of course, 0 is a function of both I and a, and hence must 

o b 

be determined separately for each computed point on these graphs). Note the 
rapid performance degradation as the rms surface deviation increases, even for 
the relatively low background level represented by = 0.04 counts/pico-sr 
(this level corresponds roughly to an average of 0.13 counts for a f ield-of- 
view half-angle of 0 = 1 yrad). 
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4.1.3 


Focal-Plane Detector Arrays 


In high background radiation environments , receiver performance can 
often be improved by using several detector elements to partition the field-of- 
view, followed by the MAP processing developed in Section 3. For the cir- 
cularly symmetric signal intensity distributions under consideration, the 
detector elements take the form of concentric rings designed to collect signal 
energy from roughly equal-intensity contours in the focal plane. The output 
of each detector element is multiplied by the appropriate weighting factor 
(Equation 3.6b) and summed to obtain the decision functions defined by 
Equation 3.5. Now, Equation 3.11 must be evaluated to determine receiver 
performance (recall that Equation 3.11 applies only to equilikely, equal-energy 
M-ary PPM signals). We have seen in Section 2 that the signal intensity 
distributions under consideration here can be roughly divided into diffracted 
and scattered components. Therefore, it seems reasonable to divide the field- 
of-view into two regions by using concentric detector elements matched to each 
component. The inner detector element (of radius 0^f) collects mainly 
diffracted signal energy, while the outer detector ring (of outer radius 
and inner radius 0^f) collects scattered signal energy. Best performance 
is obtained by finding the two radii that minimize the probability of error, 
subject to the constraint 0^ > 0^. (Of course, further improvements are 
possible by using more than two detector elements; however, for the type of 
intensity distributions under consideration, the additional improvement is 
generally not significant.) The improvement over the optimized single- 
detector configuration can be seen in Figures 4-4(a) and (b). (Note that in 
Figure 4-4(a) , M = 2.) The greatest relative improvement is obtained in 
Figure 4-4(a) , where most of the signal energy is in the scattered fields, but 
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Figure 4-4. Minimum Error Probability as a Function of K g0 for 
Optimized Single-Detector (N = 1) and Partitioned 
Field-of-View (N = 2) with: (b) High Background Level, 
M = 32 
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Figure 4-4, Minimum Error Probability as a Function of K go for 
Optimized Single-Detector (N = 1) and Partitioned 
Field-of-View (N = 2) with: (c) Fractional Signal 
Power vs. 0 for Cases Considered in Figures 4-4(a) 
and (b) 
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a significant amount of essentially noise-free signal energy remains in the 
diffracted field. A lower bound on the achievable performance is also shown 
for comparison ( = 0). Note that significant signal power savings can be 
realized by using two detector elements (N = 2) relative to the single-detector 
case (N = 1), but the noiseless performance bound cannot be reached, even as 
the number of detectors grows without bound. Less favorable situations are 
illustrated in Figure 4-4(b) , where the noise level, and the distribution of 
diffracted and scattered signal components are such that no significant 
advantage can be attained over the optimized single-detector case. The signal 
intensity distributions corresponding to the performance curves of Figures 
4-4(a) and (b) are shown in Figure 4-4(c) . 


4.2 PERFORMANCE BOUNDS AND APPROXIMATIONS 

The bounds and approximations developed in Section 3 can often be 

used to obtain quick answers without extensive numerical computation. The 

approximate expression given in Equation 3.17 yields excellent results for both 

single- and multiple-detector receivers as the background intensity becomes 

sufficiently great. This approximation is based on the fact that the envelope 

of the Poisson probability density approaches a Gaussian density as the average 

count-rate grows without bound. It is apparent from Figure 4-3 that for the 

single detector case (N = 1), the agreement between the exact solution and the 

Gaussian approximation is very good when > 10, but the approximation tends 

to become unreliable for small background counts (i.e., when 1). When 

two (or more) detectors are employed, the Gaussian approximation continues to 

yield good results as long as all detector elements observe high background 

counts. In Figure 4-6, condition "A" is defined as = 10, = 5, 

K = 2 K ., and condition "B" asfG =2, K, 0 = 1, K 0 = 2K , . Note 
s2 si* bl ’ b2 s2 si 
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Figure 4-6. Symbol Error Probability Approximations for a Partitioned 
Field-of-View with: (b) M = 32 
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that the exact result and the Gaussian approximation are in good agreement 
under condition "A" for M = 32 as well as for M = 2, but that the approximation 
tends to become erratic under the weaker background condition, "B." Thus, care 
must be exercised when using this approximation in regions where the Gaussian 
assumption may not apply. 

The behavior of the binary upper bound defined in Equation 

3.20, is shown in Figure 4-5(a) as a function of K g , for the case of inter- 
mediate (K^ =1) and large (K^ = 10) background counts. The M-ary bound 
(M - 1) is shown in Figure 4— 5 ( b ) for M = 32. It is apparent that these 
bounds are tight at low error probabilities, for intermediate as well as large 
background counts. The approximation to the binary bound defined in Equation 
3.21b uses the first term of the asymptotic expansion of I (i|;); therefore, 
it is accurate as long as >> 1 , a condition that is guaranteed to hold if 
>> 1 . (For the high-background case shown in Figure 4-3, the approximation 
cannot be distinguished from the bound in the region K g > 5.) However, the 
low background count approximation to the binary bound (Equation 3.21a) should 
only be used in situations involving extremely low background counts , as it 
tends to become inaccurate in intermediate-background cases. 
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SECTION 5 


CONCLUDING REMARKS 


The possibility of using "rough” reflectors for the purpose of receiving 
optical signals has been examined in the previous sections. It was shown that 
reflector surface quality has little effect on receiver performance in the 
absence of background radiation, as long as the receiver's f ield-of-view can 
be increased sufficiently to collect most of the scattered signal energy. 
However, background radiation may cause serious deterioration in receiver per- 
formance when rough reflectors are employed. This problem can be ameliorated 
by partitioning the receiver's f ield-of-view to obtain independent samples of 
the distorted signal field, and resorting to more complicated schemes for pro- 
cessing the field samples to improve receiver performance. It was found that 
under favorable circumstances adequate performance can often be maintained 
using nondiffraction-limited receivers, possibly leading to significant savings 
in cost and receiver complexity. Bounds and approximations to the exact error- 
probability expressions have also been developed to aid in system design 
without the need for extensive numerical computation. 
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APPENDIX A 


In this appendix, a useful form of the conditional probability 
til 

density of the i weighted count, given , is derived. Let denote 

the sum of all elemental counts weighted according to hypothesis "i": 


z . 
1 


L 


£-1 


N 

£ 

n=l 


k *n 


(A.l) 


Since the elemental counts are independent, the characteristic function of z^ 

can be expressed as the product of the characteristic functions of the 

individual weighted counts. If the counts are governed by hypothesis "q," then 

the characteristic function of z . , given H , becomes 

1 q 


Y [ exp ( J “ ”i„) - ‘1 J 

4=1 n=l ) 

The conditional probability density of z^, given H^, can be expressed as 

the inverse transform of d> (u)|H ): 

z . 4 


(a| 


V 


■if 


<t>„ (a> 


H )e^ (0a dw 

q 


(A. 3) 


Note that p (a 
i 

convolution of a finite number of discrete densities. Therefore, this condi- 
tional density can always be expressed as 


|H ) is a discrete density, since it can be viewed as the 

1 q 
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(A. 4) 


P (<x|H ) = 

Z, q 


E p (i, H ) 6 (a - o£) 
k q k 


k=0 


where p (i,H ) represents a nonzero probability mass defined over the point 
k. q 

a = a^. With no loss in generality, we can assume the ordering 


n 1 / 1 / i 

0 = <*0 < . . . 

so that for any M i," {a^} is an ordered set of points along the nonnegative 
real line. The value of the probability mass over some point depends on 
the distribution of signal and background energy over the detector array. Given 
, the values of p^(i,H^) can be computed for any k, once the distribution of 
signal and background energy is specified. 


A- 2 



APPENDIX B 


This appendix is devoted to the proof of the inequality shown in 
Equation 3.18. By symmetry, it suffices to prove that 


P (E | H ) < (M-l ) P (E | H ) 
1 M 1 — 12 1 


(B.l) 


Write the conditional error probability as 


lV E| V ' P M1 + P M2 


(B.2) 


where 


and 


P M1 = Pr {k 2 > k l or k 3 > k l * *• or ^ > k i • H P 


P M2 = Pr { comm ^- tt: *- n S an error when attempting to resolve ties 
among maximal counts including the signal count | . } 


If 

P M1 < ( m “OP 21 (B.3a) 

and 

P M2 1 (M_1)P 22 (B.3b) 

then it follows that Equation B.l is true. But Equation B.3a is true because 
the probability of the union of events is upper-bounded by the sum of the prob- 
abilities of the constituent events (union bound). To show that Equation B . 3b 
is true, let 
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x k = k! 


b e b 


and 


k -‘ V e'S 


E 


y k = 7 . i 

J-o 


and write P„„ as 
M2 


oo 

P M2-2 C<K » ,K s ;k> 

k=0 


M-l 

Z / r \ /M-l\ r-1 M-r-1 

(ttt) l x) x k y k 
r-1 


where 


<VV * ( W 


c < K b . K s i« =— — e 


Since l C ( , K g ; k ) = 2 P 22 » it suffices to prove that 


(B .4) 


M-l 

M-l\ r-1 M-r-1 _ (M-l) 

, r ) x k y k < —2~ 

r=l 



(B.5) 


But 



Therefore 


M-l 

M-l 

£(58 C-S 

i r-1 M-r-1 /M-l\\ ^ /M-2\ r-1 M-r-1 

' x k y k < (—)/ 4 \r-l/ x k y k 

r-1 

r-1 


/M-l\ , . .M-2 . M-l 

(—) (x k +y k> * — 


Hence Equation B.l is true as asserted. 
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